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ABSTRACT 



In early 2008, the young low- mass star EX Lupi, the prototype of the EXor 
class of eruptive variables, optically brightened by over five magnitudes for a 
period of 7 months. The previous time a change of such amplitude had been 
observed in EX Lup was over 50 years ago. In this Letter we present new optical 
and near-IR high resolution spectroscopy of EX Lup during the 2008 outburst. 
We investigate the physical characteristics of the outburst both soon after it 
began and some four months later, and consider the energetics and kinematics 
observed. Emission line strengths, widths, and profiles changed significantly 
between the two observations. Also, modeling of the 2.2935 /im CO overtone 
bandhead emission suggests that an inner gap in the circumstellar gas disk around 
^ \ the star may be present and it is from the inner edge of the gas disk that the 

CO overtone emission probably arises. We derive a mass accretion luminosity 
and rate during the extreme outburst of ~2±0.5 L and ~2±0.5xl0 -7 M yr _1 , 



respectively, which suggests that this outburst was indeed one of the strongest 
witnessed in EX Lup, yet not as strong as those observed in FU Orionis stars. 

Subject headings: stars: individual (EX Lupi) -- circumstellar matter -- stars: 
formation - - stars: variables: T Tauri, Herbig Ae/Be - - accretion, accretion 
disks 
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INTRODUCTION 



Eruptive outbursts in young low-mass stars are scarce events that are thought to be 
the result of increased accretion through a circumstellar disk (Hartmann & Kenyon 1996). 
EX Lupi is the prototype of a class of young stars (termed EXors) showing such outbursts 
(Herbig 1977, 1989, 2007). It has been observed to exhibit sporadic brightness increases of 
up to a few magnitudes which typically last a few weeks before fading back to a quiescent 
state. Such a 'characteristic outburst' occurred in 1993 and was studied spectroscopically by 
Lehmann, Reipurth, & Brandner (1995) and Herbig et al. (2001). In contrast to these smaller 
characteristic outbursts, EX Lup was observed in 1955 to exhibit a larger outburst resulting 
in an optical brightness increase of ~5 magnitudes. This event, which we herein term an 
'extreme outburst', was until recently unique in the photometric annals of EX Lup. However, 
in early 2008 a second such event took place. Since it is unclear whether extreme outbursts 
are initiated by the same physical mechanism as their smaller more frequent counterparts, 
the advent of the 2008 eruption presented us with the opportunity to perform a detailed ob- 
servational study of an extreme outburst. Below, we present the results of this investigation 
utilizing high spectral resolution optical (Keck I/HIRES) and near-IR (Keck II/NIRSPEC) 
observations. 



2. OBSERVATIONS AND DATA REDUCTION 

Observations of EX Lup (16:03:05.5, -40:18:25, J2000) were made on the Keck I and II 
10- meter telescopes on UT 2008 January 23 (J23), 2008 February 23 (F23), and 2008 May 
23 (M23). Two instruments were used: HIRES (Vogt et al. 1994) on Keck I and NIRSPEC 
on Keck II (McLean et al. 1995). Optical high-resolution spectra were obtained with HIRES 
on J23 and M23. On J23, the HIRESr configuration was used with the KV380 filter and a 
slit width of 0'.'861 resulting in a spectral resolution of R~48,000 (i.e. a velocity resolution of 
~7 km s _1 ) and a wavelength coverage of ~3850-8350 A. On M23, the HIRESb configuration 
was used with the WG335 filter and the same slit resulting in a similar spectral resolution 
and a wavelength coverage of ~3500-6300 A. We are greatly indebted to G.H. Herbig for 
making the J23 HIRES spectra available. Near-IR (NIR) high-resolution spectra (R~18,000, 
velocity resolution ~17 km s _1 ) were obtained with NIRSPEC on F23 and M23. 
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3. RESULTS 

3.1. Optical Light Curve 

Fig. [T] shows the optical light curves of EX Lup over the 1955 (top), 1993 (middle), and 
2008 (bottom) outbursts, taken from the AAVSO database. The 2008 outburst plot shows 
that the J23 HIRES observations were acquired soon after the extreme outburst occurred, 
some 8 days after we had been notified of the event (A. Jones, private communication to 
G.H. Herbig) and while the source was still in a rapid brightening phase with my ~10.5 or 
~3.5 magnitudes brighter than in quiescence (my ~14). Our first NIRSPEC observations 
were taken one month later when EX Lup had passed its peak brightness of my ~8 and had 
my ~9.5. Our second set of HIRES and NIRSPEC observations were taken simultaneously 
three months later (M23), during the relatively constant 'plateau phase' of the outburst 
which lasted ~170 days. At this time EX Lup had my ~10.2 and was close to the start of 
a period of instability, about two months prior to the sharp decline back to quiescence. We 
note that both the 1955 and 2008 extreme outbursts lasted approximately the same length 
of time although they are somewhat different in behavior; the 1955 eruption started slowly 
and built up to a maximum (after ~160 days) while the 2008 event appeared to have a 
more rapid rise to maximum. The light curve suggests that the 1993 event was a smaller 
characteristic outburst. 



3.2. HIRES Spectra 

In Fig. [2] we show three wavelength regions from the J23 HIRES spectrum, selected to 
cover important atomic features, specifically, H/3 (4861.3 A top), He I (5875.7 A) and Na I 
(5890.0 and 5895.9 A middle), and He I (6678.2 A) and Li I (6707.8 A bottom). All but Li I 
are seen in emission (which is seen in absorption). Na I and H/3 both possess blue-shifted 
absorption, with the emission lines being broad (full width 10% intensity ~750 km s _1 ) and 
the blue-shifted absorption removing flux well below the continuum with a minimum at ~- 
120 km s _1 1j The blue-shifted absorption on H/3 appears smooth with one broad minimum 
at -115 km s . Similarly, the blue-shifted absorption on the Na I D lines have minima at 
-129 and -122 km s _1 , respectively. Fe I, Fell, and He I all exhibit composite profiles made 
up of both broad and narrow emission components (BC and NC, respectively; Hamann & 
Persson 1992). The NC and BC peaks lie within ~10 km s _1 of their rest wavelength. 



1 All radial velocities are relative to the EX Lup rest velocity. EX Lup has Vheiio — +23, +28, and 
-2 km s -1 , on J23, F23, and M23, respectively. 
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Four months after the first HIRES observation, on M23 (see Fig. [3]), much of the same 
spectral structure was still visible although some changes appeared to have occurredj. The 
He I line still possessed NC and BC features although the radial velocity offset of the NC 
from the line rest wavelength had shifted to the blue. Also, the H/3 profile had changed and 
exhibited a double peak at -77 and +100 km s" 1 with an absorption minimum at +42 km s _1 . 
However, the width of the H/3 emission was similar to that found four months previously. 
The profiles of the Na I D lines were also quite different and became somewhat more com- 
plex. Both lines remained in emission and possessed blue-shifted absorption components 
although the latter were much weaker than in the J23 spectra. Additionally, the blue-shifted 
absorptions showed radial velocity structure (see Fig. [3]) extending to -277 km s _1 on M23 
and -212 km s _1 on J23. Finally, superimposed on both of the broad Na I emission lines 
were narrow absorption features with a radial velocity offset of -4 km s _1 . 

We have performed a best fit to the NC and BC structures present in He I and Fe II 
using a two-component gaussian emission profile model, iterating over three parameters for 
both gaussians, specifically, the wavelength of peak emission (A pea £;), the peak flux (F pea k), 
and the gaussian FWHM. The results of this procedure are shown in Fig. HJ Between the 
two dates we find that, i) for both lines \ ea k of the NC and BC gaussians remained within 
~10 km s _1 of the line rest wavelength, ii) the He I NC FHWM did not change significantly, 
remaining at ~36 km s _1 , while the Fe II NC FWHM became smaller by a factor of 2 (14 and 
7 km s _1 ) and at both times was significantly smaller than the He I NC FWHM, in) the He I 
BC FWHM declined (from 220 to 180 km s" 1 ), while the Fe II BC FWHM remained around 
the same (at ~170 km s _1 ), iv) the He I ~F pea k of both the NC and BC increased slightly (by 
~20%). The He I equivalent width (EW) on the two dates was EW(He I) = -3.5 A (J23: 
EWtvc = -0.69 A, EWbc = -2.80 A, M23: EWtvc = -0.82 A, EW BC = -2.70 A). Similarly, 
EW(Fe II) on J23 was -5.8 A (EWtvc = -0.21 A, EW BC = -5.60 A) and on M23 was -7.9 A 
(EWtvc = -0.10 A, EW BC = -7.80 A). 



3.3. NIRSPEC Spectra 

Data for three spectral regions taken from the NIRSPEC observations are shown in 
Fig. [5j Here, results from both observing dates (F23 and M23) are plotted for the regions 
covering Br7 at 2.166 fim (top-left) and the v=2-0 CO overtone bandhead at 2.2935 /iin 
(top-right). The bottom-left plot covers the Na I doublet lines at 2.206 and 2.208 /im from 



2 In the M23 observations, Li I (6707.8 A) and He I (6678.2 A) did not fall on the detector. 
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F23 onljo The Br7 line appears very similar in J23 and M23 with a EW(Br7) of -17 A. It is 
broad (FWHM~220 km s" 1 on F23, 190 km s" 1 on M23), and is slightly red-shifted with the 
Xpeak emission offset by +5 and +32 km s _1 on F23 and M23, respectively. The Na I doublet 
lines are also broad with FWHM of 177 and 155 km s _1 on the one date they were observed. 
The most striking change occurs in the spectral region containing the 2.2935 /iin v=2-0 CO 
overtone bandhead. In the data from F23, the bandhead is strongly in emission with a 
EW(CO) = -21 Ac| and exhibits the characteristic signature of being significantly velocity 
broadened (i.e. a 'blue-hump' is present, see for example WL 16 in Carr et al. 1993). In 
the data from M23, however, the emission is only barely seen (EW(CO) = -4 A) although 
velocity broadening of the bandhead still seems to be present. In the bottom-right plot of 
Fig. [5] we have modeled the CO bandhead emission profile using the analysis from Dent 
& Geballe (1991) and Carr & Tokunaga (1992) previously employed in Aspin, Reipurth, 
& Lehmann (1994) and Aspin & Greene (2007). The details of the model are presented 
in the aforementioned papers and, for brevity, here we merely present our findings. The 
CO bandhead emission is best fitted using a velocity profile originating in a circumstellar 
disk. This produces a good fit to the blue-hump and the overall shape of the bandhead. 
We investigated a wide range of the model free parameters, specifically stellar mass (M*) 
over the range 0.1-2.0 M , the emission region temperature (Tco) from 1500 to 6000 K, the 
optical depth of CO {tco) from 0.01 to 5.0, and the inner (r in ) and outer (r out ) gas disk radii 
between l.OxlO 9 and l.OxlO 11 m (i.e. ~1.4-140 R* or 0.007-0.7 AU). The best fit to the 
spectrum gave values of M* = 0.6 M , T co = 2500 K, t co <0.1, r in = 1.2 xlO 10 m (~17 R* 
or 0.08 AU), and r out = 2.0 xlO 10 m (~29 R* or 0.13 AU). Typical associated uncertainties 
on these values should be considered to be ~20%. For this value of r in and with R* = 1.6 R 
(as estimated by Gras- Velazquez & Ray 2005 and Sipos et al. 2010), gas at 2500 K would 
require a heating source on the stellar surface of ~8000 K similar to what was found for 
surface accretion shocks in, for example, V1647 Ori (Briceno et al. 2004). The Keplerian 
velocity at 0.08 AU for M* = 0.6 M Q is ~80 km s _1 . Due to the complex relationship between 
the model free parameters, the uniqueness of the above best-fit model is unclear. However, 
we believe that our investigation of the aforementioned free-parameter space results in values 
for the above physical properties that are reasonable for this star. 



3 The Na I spectral region was not included in the M23 dataset. 
Calculated from 2.2925 to 2.2975 /j,m. 
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4. Discussion 

4.1. Photometric variability 

The 2008 extreme outburst of EX Lup appears to be unlike any other observed in this 
source since 1955. Both extreme outbursts exhibited very large optical brightness increases 
(5-6 magnitudes) and, although information about the earlier event is scarce they appear, at 
least superficially, to be quite similar. Between 1955 and 2008, EX Lup has undergone nu- 
merous smaller characteristic outbursts with Amy typically 1-3 magnitudes. The durations 
of the extreme and characteristic events appear very different; the former last 200+ days 
whereas the latter last perhaps only a few weeks (see the AAVSO archive for more details). 
Photometrically, therefore, there seems to be an inherent difference between extreme and 
characteristic outbursts. 



4.2. Spectroscopic Variability 

The HIRES spectra of EX Lup taken between 1998 and 2004 (presented in Herbig 2007) 
and covering both the quiescent state and the characteristic outburst phases, showed narrow 
emission features together with an underlying atomic and molecular absorption spectrum 
typical of an ~M0 dwarf star. During the 2008 extreme outburst, however, the HIRES 
spectra appeared markedly different with no obvious sign of an underlying photosphere. 
Comparing our data taken during the 2008 extreme outburst (shown in Figs. [2] and [3]) with 
previous observations we find that whereas Herbig et al. (2001) found strong emission lines 
of H I, He I, He II, and Ca II with no evidence for blue-shifted absorption nor Fe I and Fe II 
lines, we observe both. Also, the characteristic outburst phase spectrum from Lehmann, 
Reipurth & Brandner (1995) possessed H I, He I, He II, Ca II, lines together with weak Fe II 
lines, exhibiting inverse P Cygni (IPC) profiles (red-shifted absorption) on H I, Ca II and 
perhaps even Fe II which disappeared around 7 weeks later. Herbig (1950) noted that IPC 
profiles were also present in his 1949/1950 spectrograms taken several years before the earlier 
extreme outburst. Lehmann et al. (1995) interpreted these IPC profiles as originating in 
cool infalling gas. Our HIRES spectra, which shows blue-shifted absorption, suggest that 
we are observing the opposite, i.e. cool outflowing gas. We speculate that perhaps it is only 
during extreme outbursts that strong and dense enough stellar/disk winds are generated to 
create the blue-shifted absorption component. 

Two-component (BC and NC) emission line profiles, similar to the ones seen in our 2008 
EX Lup spectra, were previously reported on the stronger Fe II and Ti II lines in EX Lup in 
data taken in June 1998 (Herbig 2007). Another young T Tauri star to show such structure 
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is DR Tau (Beristain, Edwards, & Kwan 1998). The interpretation suggested for this source 
was that the NC of the lines arose in a compact region close to the stellar surface (and 
hence not significantly rotationally broadened) as in a post-shock region at the base of an 
accretion flow. The BC was considered to originate in a region with considerable radial 
velocity structure such as in a magnetospheric accretion funnel or in the hot inner region of 
the circumstellar disk. The different excitation requirements for He I (E upper .=23 eV) and 
Fe II (E Mpper =5-6 eV) suggests that these lines are created in different physical locations, 
making the similarity of structure somewhat peculiar. Any model proposed to explain the 
formation of these lines and the creation of the NC+BC structure would need to explain this 
energy discrepancy together with the temporal changes observed between J23 and M23 i.e. 
i) a decline in the strength (i.e. EW) of the Fe II NC by a factor 2 coupled with a ~40% 
increase in the BC strength, it) a ~20% increase in the strength of the He I NC coupled 
with an unchanged BC strength, and Hi) the changes in kinematic properties (e.g. FWHM, 
radial velocity offset) shown in Fig. HI 

The changes in the blue-shifted absorption profiles between J23 and M23 suggest that 
the wind from EX Lup evolved between the two observations. The blue-shifted absorption 
features present on J23 and seen in the H, Na, and Ca lines, varied somewhat differently for 
each element. Between J23 and M23, i) the blue-shifted H/3 absorption moved redward and, 
in fact, became somewhat red-shifted, ii) the blue-shifted Ca I absorption disappeared, and 
Hi) the deep and broad single component blue-shifted absorption on the Na I lines became 
shallower and more complex in structure, and extended to much higher radial velocities. A 
sophisticated radiative transfer code is required to model such variations. 

The changes in the CO overtone bandhead observed between J23 and M23, together 
with the modeling presented above, suggest that the hot gas in the inner regions of the 
accretion disk either cooled significantly between the two dates, or was mostly accreted onto 
the star, or the UV photoexcitation conditions changed (see Brittain, Najita, & Carr 2009). 
The unchanged Br7 flux over this period implies strong on-going accretion. 



4.3. Accretion Luminosity and Rate 

On J23, F23, and M23, EX Lup had my values of 10.5, 9.5, and 10.2, respectively. 
Juhasz et al. (2010, in preparation) measured the K-band magnitude of EX Lup to be 
Hift-=7.5 on UT 2008 April 20. On that date EX Lup had a V magnitude of ~9.5 and hence 
a V-K=2.0. If we assume this V-K value for F23 and M23 then on those dates m^ was 
7.5 and 8.2, respectively (a mean of 7.85). The Br7 emission line fluxes therefore decreased 
by 2x between J23 and M23 since on both dates the EW of Br7 was -17 A. If the optical 



continuum emission and the hydrogen line fluxes are indicative of the mass accretion rate 
(e.g. Gullbring et al. 1998; Muzerolle, Calvet, & Hartmann 1998; Muzerolle, Hartmann, & 
Calvet 1998) then the accretion rate changed by only a factor 2 over this four month period. 

We can convert line and continuum luminosities to accretion luminosity (L acc ) to esti- 
mate the mass accretion rate (M). Following Aspin, Beck, & Reipurth (2008) and adopting 
m^-=7.85, a flux corresponding to zeroth magnitude at 2.2 jivcv of 4.35 xlO -10 W m -2 yum -1 
(Tokunaga & Vacca 2005), a visual extinction Ay ~0, r in =17 R* (see above), L*=0.5 L 
(derived from T*=3800 K and R„,=1.6 R , again see above), and a distance of 140±20 pc 
(Hughes et al. 1993), we obtain values of L acc 1.6±0.6 L Q and M=1.5±0.5xl0~ 7 M yr -1 . 
Following Valenti et al. (1993) and Gullbring et al. (1998), the change in optical contin- 
uum emission in extreme outburst (AV=4.5 magnitudes) is converted into L acc ~2.3 L and 
M ~2.1xl0 -7 M Q yr -1 . We therefore adopt extreme outburst values of L acc =2±0.5 L Q and 
M=2±0.5xl0~ 7 M Q yr -1 . Sipos et al. (2009) estimated the quiescent accretion luminosity 
and rate of EX Lup to be L acc ~0.0037 L and M ~4.0x 10~ 10 M Q yr -1 using Pa,# line fluxes. 
Since EX Lup is veiled even in quiescencqj we adopt veiling-corrected values for quiescent 
L acc and M of 0.05±0.02 and 6±3xl0~ 9 , respectively. 

The accretion luminosity and rate during the 2008 extreme outburst of EX Lup are 
therefore ~40x larger than their quiescent values. However, these increases are still ~50x 
smaller than those found in FU Orionis eruptions (Hartmann & Kenyon 1996). 

In conclusion, the 2008 extreme outburst of EX Lup possessed a number of features 
that distinguished it from the more common characteristic outbursts and appears to have 
had more in common with the previous extreme outburst that occurred in 1955. 
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Fig. 1. — Optical light curves during the 1955 (top), 1993 (middle), and 2008 (bottom) 
outbursts. All plots cover 400 days starting from shortly before the outbursts began. The 
mean 'quiescent' optical brightness of EX Lup is indicated by the horizontal dashed lines. 
All data points are taken from the AAVSO database. The vertical dot-dashed lines in the 
bottom plot indicate the times when the HIRES and NIRSPEC spectroscopic observations 
were acquired. The dates shown in the middle plot are those when the Lehmann et al. (1995) 
optical spectra were taken. 
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Fig. 2. — Optical HIRES spectra taken on UT 2008 January 23 close to the start of the 
outburst. Three wavelength regions are shown containing H/3 (top), He I and Na I (middle), 
and He I and Li I (bottom). Heliocentric radial velocities are given for the major features of 
the atomic lines. Equivalent widths (W) of the important features are shown. 
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Fig. 4. — HIRES spectra of the Fe II line at 4233.2 A (top, grey lines) and the He I line at 
5876.1 A (bottom, grey lines) overlaid with model gaussian profiles. The best fit parameters 
are given in the legends. The NC gaussian is shown as dot-dashed lines while the BC gaussian 
is shown as dashed lines. The total model profile is the solid (black) lines and the residuals 
between the observed and model profiles are shown by the dotted lines. The heliocentric 
rest wavelength of the line is shown as a vertical solid line. 
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Fig. 5. — NIRSPEC spectra taken during the extreme outburst of 2008. In two of the four 
plots, spectra taken at two different epochs are shown offset slightly for clarity (UT 2008 
February 23 and May 23). Top-left is the Br7 emission line at 2.166 /mi. Bottom-left plot 
is the Na I emission lines at 2.206 and 2.209 /tm. Top-right plot is the v=2-0 CO overtone 
bandhead at 2.294 /im. Bottom-right is the CO bandhead profile from UT 2008 February 
23 (black solid line) together with the best-fit disk emission model (grey solid line, offset 
for clarity). The residuals between the observations and model are also shown (dotted line). 
The best-fit model parameters are described in the text. 



